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The dynamical scaling properties of selfavoiding polymerized membranes with internal dimension 
D are studied using model A dynamics. It is shown that the theory is renormalizable to all orders 
in perturbation theory and that the dynamical scaling exponent z is given by 2 = 2 + D jv* . This 
result applies especially to membranes (73 = 2) but also to polymers (J) = 1). 



Polymers and polymerized flexible membranes show in- 
teresting statistical properties. Polymers have been in- 
vestigated since long time, both static and dynamic [0. 
Selfavoiding polymerized membranes have attracted a 
remarkable interest during the last years. Their static 
properties have been studied numerically [0-p| experi- 
mentally [^pj and analytically [p-13|. Dynamics has 
been regarded using scaling arguments for polymers 
|14| , |l5| and membranes [|l6| . For polymers, a renormal- 
ization group analysis has been performed at 1- 17-22|] 
and 2-loop |2^] order. 

For membranes, the analytical approach, inspired from 
polymer theory [|l|, relies on renormalization group and 
£-expansion methods. It was initiated in [^ 10|, where 
it was used to perform calculations at 1-loop order. Its 
consistency to all orders in perturbation theory has been 
established in |ljj2J]. Recently, 2-loop calculations have 
been performed, which give reliable results for all embed- 
ding dimensions |l^ , p^ . 

In this letter we address the question of the dynam- 
ics of such membranes, always including polymers as a 
special case. 

The membrane is modeled by a continuum model a 
la Edwards: the embedding of the Z3-dimensional mem- 
brane in d-dimensional bulk space is described by the 
mapping x G IR — > r(x) € IR . The Hamiltonian, 
which describes the static properties of the membrane. 



n[r]^l f{Vrix))'+b f fs'{r{x)-r{y)) , (1) 

<y X J xJ y 

where b is the coupling constant. The dynamics of the 
membrane is given by the Langevin-equation 



f{x, t) = ~X + C{x, t) 

dr[x,t) 



(2) 



which models the purely diffusive motion of the mem- 
brane (model A in the terminology of p^]). The Gaus- 
sian noise Ci^, t) has the correlation 



C(a;, t)({x',t') = 2X5" {x - x')5[t - t') . 



(3) 



This Langevin equation can be formulated using an ef- 
fective field theory p6|] with action 



J [r, r]^ J^ Jj{x, t) l^r{x, t) + A^^^ 



\f{x,tf 



(4) 



Expectation values are calculated by integrating over all 
fields r and f. Prepoint-discretization is used and the 
average over the noise has been taken. (This generated 
the term quadratic in f.) 

Perturbation theory is performed by expanding about 
the Gaussian theory. We use the free propagator G and 
correlator C in position-space 



C{x,t) 



1 /I 



^^{r{x,t)-r{0,Q)f 

|xpT-i(-i^)x 
/ 1 /4A|t|V 



(5) 



v \ x'^ 



ds 



-^(e 



1 



G(x,<):=-(r(x,t)r~(0,0))o 

= e(i) (47rA|t|)-^/2 g-xV4A|t| ^^^2 - D) , (6) 

where Sd is the volume of the unit sphere in D dimen- 
sions, 

Sd = ^^, (7) 



and 



T{D/2) 



2-D 



(8) 



Some normalization-factors have been absorbed into the 
measure in order to have 



C{x,t)Ki \x\ 



2-D 



for x^ > A|t| 



(9) 



(For details compare appendix A of jl^.) If x^ is much 
smaller than A|t|, the correlator approaches the finite 
value 



c(-') = iwi+°<^^'' 



(10) 



We should mention that the propagator is simply related 
to the time-derivative of the correlator 



G{x,t) = e{t)^C{x,t) 
A 



fii) 



We are now in a position to construct the perturbation- 
theory. The interaction vertex is: 

w- . := 2 / f{x, t){ik) e'Hr{x,t)-r{y,t)) (^2) 

Jk 

The perturbative expansion of an observable O can be 
written as 



(^)-NormE(^/(C.^ 



(13) 



where the normalization Norm has to be chosen so that 
(1) = 1 and the integral is taken over all arguments of 
the interaction vertex. We claim that divergences only 
occur at short distances and short times. To prove this 
look at a typical expectation value 

(^(y mM • )q ^ / ^ / Ja\Xi XjYn H *mj "^i , k,Ya) X 

y^e-E.,,Q^^>''''^ ^ (14) 



The exponential in (|lj) now is 



d^p f du 



^ hk,Q., = (2 - D)Sd J -^ J 



where each contribution consists of a function /„, which 
is a product of propagators, correlators and fc's and an 
exponential factor, with 



Q^ 



O yXi Xj , Li Lj j . 



(15) 



fa is a regular function of the distances. Divergences at 
finite distances can only occur if Qij is not a positive 
form. We will show that Qij is a positive form for all ki 
which satisfy the constraint 



Y^h^o 



(16) 



This constraint always holds, see equation (g^. For equal 
times it is just the statement that the Coulomb energy 
of a globally neutral assembly of charges is positive. One 
simply identifies C with the Coulomb-propagator and ki 
with the charges. In the dynamic case, write 



*-<--)^"/i^/£ 



2A 



uj^ + [Xp^f 



(f^ip(xi~Xj)+iui(ti-tj) _ -^ /^ij\ 



2A 



{2n)D J 2^ ^2 



d^p f du 



'^-">*-/(^/ 



(27r)^ J 2tt 



UJ^ 
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^fc./cj L^pi^^-^,)+Mu-tj) _ A 



j:^ 



^ipxi-\-ujti 



(18) 



To get the second line, equation ( [iq ) has been used. Note 
that again due to equation (|l6|) , the integral is ultraviolet 
convergent and thus positive. It vanishes if and only if 
the charge-density, regarded as a function of space and 
time, vanishes. This is possible if and only if endpoints of 
the dipoles (which form the interaction) are at the same 
point in space and time. No divergence occurs at finite 
distances. To renormalize the theory, only short distance 
divergences have to be removed by adding appropriate 
counter-terms. In addition, as the divergences occur at 
short distances, they can be analyzed via a multilocal 
operator product expansion (MOPE) |ll|. For a detailed 
discussion of the MOPE and examples see |jri|,||,|l|,0 . 

We show now that the counter-terms which render the 
static theory (m finite are also sufficient for the dynam- 
ical case (H). As an illustration we first calculate the 
1-loop counter-terms for the renormalization of the field 
and of the coupling-constant. The first singular configu- 
ration appears, when both ends of the interaction vertex 
(O) are contracted towards a single point. The leading 



term (MOPE-coefficient) of this expansion is (the dotted 
line indicates points which are contracted): 



(19) 



= 2 / -[r{x,t){ik)]e- 

Jk 

x[{ik){r{x,t)~r{y,t))] + 
We now Taylor-expand r{y,t) as 
T{y, t) = r{x, t) + {y~ x)Wr{x, t) 

+ l[{y-x)S/fr{x,t) + 0{\x~y\'') (20) 

The leading term in equation (n9[) is 

r[x,t)(~Ar(x,t)) ;; ;— e 



D d 



f(x,t)(-A)r(x,t) I — ) ix-y 



D-ud 



(21) 



Denoting with 



wf the MOPE coefficient of 



proportional to -^ 
be written in the form 



Vu\i— 



f{x,t){—A)r{x,t), this can 



Q 



+ 



where 



^l=-2i^l^ 



\D-ud 



(22) 



(23) 



is the static MOPE-coefficient [^JJlJ. This imphes that 
the counter-term for the wave-function renormahzation 
is the same as in the static case. Let us now regard 
the counter-term for the couphng-constant renormahza- 
tion. Using the techniques of pl|-p^, we obtain for the 
contraction of two interaction vertices to one interaction 
vertex: 



-{G{x,t)+G{y,t))x 



x[Cix,t)+Ciy,t)r''/'-' (24) 

X and y are the distances of the contracted endpoints of 
the dipoles, t is their time-difference. The trick is now to 
write this expression with the help of equation (O) as 



±e(t>0)^^[Cix,t) + Ciy,t)] 



-d/2 



(25) 



In performing the perturbation-theory, we have to inte- 
grate over all times. If we use, as is usually done, no 
cut-off in the time-direction, the time integral will sim- 
ply give the value of the function at its lower bound: 



dt 




(26) 



(The r.h.s. is the counter-term of the static theory, see 
[^-|l3[.) We easily convince ourselves that this relation 
implies the same countcrterm as in the static case, if we 
take care of the additional combinatorial factor 2 for the 
time-ordering of the interaction vertices. 

One knows from general arguments that the diver- 
gences associated to short distances in space are removed 
by the static counter-terms |2^ . We now use the fact that 
the static theory is renormalizable |jll| . This implies that 
new divergences can only appear for short times. We 
therefore have to analyse all possible divergences of this 
type. Using the MOPE, the most relevant divergences 
are associated to the operator 

■^ ^f{x,t)f{x,t) . (27) 

We now regard a general contraction of n dipoles towards 



(28) 



In order to obtain the operator ™^ , one has to contract 
all fields r and f except of the field f (z, t) with the largest 
time-argument (all other contractions give 0). One also 
has to leave uncontracted one arbitrarily chosen field r. 
Due to the structure of the interaction (12|), the field r 
always appears in the form r{x, t — t) — r{y, t — t). So 
the contraction yields 

f{z, t) [r(x, t-r)- r{y, t - r)] M(distances) , (29) 

where M denotes the MOPE-coefficient which depends 
on the distances in space and time. Now, r{x,t — t) — 
r{y,t — t) has to be expanded about {z,t). The leading 
term has at least one spatial gradient. No term of the 
form of ( p7| ) can be constructed. Therefore, there is no 
singular contribution of this type in any order in pertur- 
bation theory and no renormahzation of ™»- is needed. 
The last at e = marginal operator is 



'{x,tf 



(30) 



Its renormahzation is given by the fluctuation-dissipation 
theorem, see below. 

We now introduce renormalized quantities according 
to 

r — y Zrn 

r= ^fn 

^ = Z\\ji 

b - bRZbZ''/^fi' (31) 

s is the dimensional regularization parameter, defined by 

e^2D -vd . (32) 

The fiuctuation-dissipation theorem states that 

e(t)|^i {r{x, t) - r(0, 0))^ ^ = A^r(x, i)f (0, 0)^ (33) 

This relation for the full expectation- values holds as well 
for renormalized as for bare quantities. We therefrom 
deduce that 



(34) 
has not to be 

(35) 



Z 

In addition, our findings that the term 
renormalized indicate that 



It follows that 



Z 



(36) 



As usual, the renormahzation group /3-function f3{bfi) 
and the full scaling dimension vibji) of r are obtained 
from the variation of the coupling constant and the field 



with respect to the renormahzation scale /i, keeping the 
bare coupUngs fixed. They are written in terms of Z and 
Z\) as 






-ebf 



l + bn^\nZ, 
2- D 






\nZ 



\^^'-^i-,'-' 



(37) 
(38) 



The /3- function has an IR- fixed point for 6|j. > 0. The 
function v[hii) is for large scales therefore given by 



u{h*) 



(39) 



We deduce that the correlation-function scales for equal 
times as 



((r(x,0)-r(0,0))2)^|x 



2v* 



(40) 



If we define the exponent z for the auto-correlation- 
function as 



((r(0,i)-r(0,0))2)^|tp/- 
the exponents z and v* are dependent: 



V* 



(41) 



(42) 



For polymers (D=l), this relation was given using scal- 
ing arguments in p4| , for membranes {D = 2) in p6[ . 
This result was followed by perturbative calculations for 
polymers in 1-loop |17|-|2^ and 2-loop [g3| order. 

It is interesting to note that for polymers and mem- 
branes (E2h can be written in the form 



y 



(43) 



where df is the fractal dimension of the membrane or the 
polymer. 

It seems attractive to use (|4^ ) in Monte Carlo simula- 
tions to improve the determination of df for selfavoiding 
membranes. This might be advantagous as z can be mea- 
sured for any point and any time, i.e. on a large statistical 
ensemble. 

In conclusion: We have shown that the purely diffusive 
motion of polymers and polymerized membranes is given 
to all orders in perturbation theory by (^3|). Experimen- 
tally however, an additional hydrodynamic interaction is 
present. This question will be addressed in a subsequent 
publication. 
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